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ABSTRACT 

The formation rate of quark matter droplets in neutron stars is calculated from 
a combination of bubble formation rates in cold degenerate and high tempera- 
ture matter. Nuclear matter calculations of the viscosity and thermal conductivity 
are applied. Results show that droplets form only in the core of neutron stars 
shortly after supernova collapse, where pressures and temperatures are high, and 
' for sufficiently small interface tension between nuclear and quark matter. Coulomb 

energies hinder formation of large droplets whereas the presence of strange hadrons 
' in nuclear matter increase the droplet formation rate. 
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^ . 1. Introduction 

■ Quark matter in neutron star has recently been found to be able to coexist 

with nuclear matter and have a rich structure in this mixed phaseS'i. This is in 
O ■ contrast to standard mo delJ where neutron stars have a pure quark matter core 

^ ' with a mantle of nuclear matter outside. The quark and nuclear matter mixed phase 



appear at much lower densities. A lattice of quark matter droplets appear already 
around twice nuclear matter saturation densities which increases the probability for 
the existence of quark matter in neutron stars. 

We shall here discuss the formation of quark matter droplets and calculate the 
^ ■ rate by which the first quark matter droplets form at the high densities in centers 

of neutron stars. (If quark matter is the true ground state of hadronic matter in 
vacuums, such strangelets would probably reside inside the stars as relics from the 



ix. 



X 



\-{ • early universe and would convert the whole star to a strange star.) We combine 



bubble formation rates of both cold degenerate and hot systems and apply transport 
coefficients of nuclear matter. To estimate the time scales for the droplet formation 
the pressure and temperatures are calculated from standard equation of states used 
in supernova core collapse for nuclear matter and a simple Bag model for quark 
matter. 

2. Droplet Free Energy 

In the core of a supernova collapse or the later neutron star the nuclear matter 
may be dense enough that it is energetically favorable to form a quark matter 
droplet. The work forming a droplet of radius R is 

Ayr 

W = R^AP + 47iaR'^ + NgAfi + Ec, (1) 
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where AP = Pqm — Pnm is the pressure difference, a the surface tension, A/i = 
fJ'QM — fJ^NM the difference in chemical potentials between nuclear and quark matter, 
Ng = ^R^uqm the number of quarks in the droplet and Ec the Coulomb energy of 
the droplet. 

It is usually assumed that the embedded droplet is in chemical equilibrium with 
its surroundings at nucleationH, i.e., 

yUn = /i« + , fip = + ^^d ■ (2) 

Olesen and Madseni consider the instantaneous collapse of the nuclear matter into 
a quark matter containing the original flavor composition. The resulting droplet 
will be out of chemical equilibrium and diffusion of various particles will take place 
until chemical equilibrium is attained. We shall take the standard approach that the 
quark matter is in chemical equilibrium with its surroundings since we expect the 
diffusion of nucleons in and out of the droplet to be at least as fast as the formation of 
a large droplet. Furthermore, we shall assume the presence of sufficient strangeness 
in nuclear matter from strange hadrons like S~, K~ , A... . With sufficient rapid 
diffusion of strangenes during droplet formation we can then also assume Hd = fis 
inside the quark matter droplet. Otherwise, only u, d quark matter droplets can 
form since strangeness production only occur slowly on weak decay time scales. 
Piire u,d- quark matter will have higher free energies and are thus harder to form. 

Let us first discuss small droplets where Coulomb effects can be neglected. Since 
the maximal temperatures reached in supernovae are much less than the chemical 
potentials T -C /i we can also ignore thermal energies in (|l|). The nucleation rate 
contains the usual Boltzmann factor in which the work to form a critical bubble 
enters. The critical bubble is determined by the maximum of Eq. (|l]) which gives 
a critical bubble radius. 



^c = |^, (3) 



at which the work is 
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The pressure difference across the boundary of the critical bubble is balanced by the 
surface tension, Pqm = Pnm + 2a /Rc. The temperatures and chemical potentials 
are the same in the two phases and therefore the critical droplet is in complete 
phase equilibrium with the nuclear matter when created. The surface tension is 
unfortunately a poorly determined parameter and we refer to Ref.i for a discussion. 
It is probably within the range a ~ 10 — 200MeV/fm^. 

The Coulomb energy, (3/5)Z^e^/-R, grows rapidly with size for a droplet with 
constant charge densitya Screening or charge rearrangement reduces the Coulomb 
energy significantly but still the critical radius and work are found to increase due 
to Coulomb effects. A more detailed analysis shows that the formation of droplets 
more than a few fm in size is severely hindered. 



3. Formation Rate 



The droplet formation rate was calculated by Langer & Turskil 



/=fl^oe-^=/^. (5) 



The "statistical" prefactor 



2 (a\^/^(R; 
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measures the phase-space volume of the saddle point around Rc that the droplet has 
to pass on its way to the lower energy state. Here C,q is the quark correlation length 
which also estimates the surface thickness. We shall assume that it is of order the 
QCD Debye screening length ~ l/qo where = {'2/n)asJ2q fJ'l = {'^/'^)c(sNgfi'^. 
Here, is the QCD fine-structure constant and Ng the number of quark flavors 
present. With ~ 0.3, Ng = 3 and p, ~ 400 MeV we have = I/qd = 0.7 fm. Note 
that Eq. (P) is only valid when -Rc~^g which is the case for reasonable estimates of 
a and AP. The "dynamical" prefactor determines the droplet growth rate and is 

^ = l^J^T + 2{^v + 0] , (7) 
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where Aw = wqm — wnm is the enthalpy difference and A, rj and C are the thermal 
conductivity, shear and bulk viscosities respectively. The dynamical prefactor (|^) 
incorporates both the low temperature result of Langer & Turski, where the viscous 
term is ignored, and the high temperature result of Csernai & KapustatI, where the 
thermal conductivity is negligible. Since both calculations were performed in the 
hydrodynamic limit linearized in the transport coefficients, the two results simply 
add upa. 

The shear and bulk viscosities and the thermal conductivity have been estimated 
by Danielewicz0 within the Boltzmann equation with an effective nucleon-nucleon 
differential scattering cross section. The contributions to the dynamical prefactor 
from the bulk viscosity and the thermal conductivity are negligible as compared to 
that from the viscosity which typically is ?7 ~ SOMeV/fm^. The result of Langer and 
Turskil include only the thermal conduction term and is therefore not applicable 
to the case of degenerate quark matter droplet formation. For the formation of 
hadronic bubbles in quark matter one can use the perturbative QCD calcuations of 
transport coefficientsEJ. Here the viscous damping dominates as well and the shear 
viscosity is even larger than in nuclear matter. 

Still neglecting Coulomb energies we obtain from Eqs. 



assuming ag = 0.3 and A^^^ = 3. Here we have written all the quantities in typical 
units Bli^: the enthalpy and pressure differences in units of 10 MeV/fm^, the quark 
chemical potential in 400 MeV, the temperature in 10 MeV, the surface tension in 
units of 20 MeV/fm^ and the viscosity in 50 MeV/fm^. The resulting rate is very 
sensitive to a, AP and T . It is rather insensitive to the prefactor and therefore also 
the viscosity. Using the prefactor T'^ in stead reduces the rate by approximately 
11 orders of magnitude which does, however, not have much effect in (|^). Here 
the huge exponent 185 arise mainly from conversion of hadronic scales (fm~^) to 
neutron star dimensions and cooling times (km~^s~^). 

The pressure difference depends on the equation of state. It increases with den- 
sity and so droplet formation is most likely in the center of the neutron star. The 
maximum temperature achieved in supernova explosions is largest in the core with 
a value around 10 MeV (depending on the equation of state and details of the su- 
pernova collapse) and it cools down to a few MeV in 10-20 seconds^!. Consequently, 
the droplets predominantly form in the neutron star center within few seconds after 

PI II rj 

the core collapse when temperatures are highest as already pointed out in (BESEj) 
and only if the surface tension is sufficiently low. 

The total number of droplets formed is the integrated rate over volume of the 
neutron star and time after the supernova core collapse 

fRNS „ poo 

N= A-Kr^dr dt I{AP{r),T{t)) . (9) 

Jo JtsN 

The rate depends mainly on the spatial coordinate through the pressure difference 
and on time through the temperature. The pressure and temperature depends 
sensitively on the equation of state of nuclear and quark matter and on the details 
of the supernova collapse. 

To convert the core of the neutron star into quark matter, whether it is in a 
mixed phase with nuclear matter or not, we must form at least one droplet, i.e. 
> 1. According to and (^) this requires 

a < 24 MeV ■ fm^^ ■ API{1 T^g , (10) 

where Pc and Tc are now the core values. This condition is almost insensitive to the 
size and cooling time of the hot core. 



4. Summary 

A formation rate for quark matter droplets has been derived applying transport 
coefficients from nuclear matter calculations. The resulting rate is very sensitive 
to the surface tension at the nuclear/quark matter interface and to the pressure 
difference and temperature. For large droplets it is important to include Coulomb 
enegies which reduce formation rate for droplets larger than a few ~ fm drastically. 
The presence of strange hadrons and sufficiently fast strange quark diffusion allow 
the formation of strange quark matter droplets which requires less work than purely 



u,d-quark droplets and thus improve the production rate. Quark matter droplets 
form only in the core of neutron stars, where the pressure difference is highest, 
and shortly after the core collapse in supernovae explosions, where temperatures 
are highest. Eq. (|10|) gives the necessary conditions for creating quark matter in 
neutron stars. The pressure difference depends strongly on the equation of state 
used for both nuclear and quark matter and with the present uncertainties in both 
and in the surface tension, the droplet formation rate cannot be reliably estimated. 
If neutron stars burn into quark starsEJ it would happen right after the supernova 
explosion, and so this mechanism is unlikely to be the cause of gamma ray bursters. 
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